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Abstract. We present a model for a corona above a ge- 
ometrically thin standard disk around a black hole. This 
corona is fed by matter from the thin disk which evap- 
orates from the cool layers underneath. An equilibrium 
establishes between the cool accretion stream and the 
hot flow. In the inner region evaporation becomes so ef- 
ficient that at low accretion rates all matter flows via 
the corona and proceeds towards the black hole as an 
advection-dominated accretion flow (ADAF). We inves- 
tigate this transition for accretion disks around stellar 
black holes. We derive a set of four ordinary differential 
equations which describe the relevant processes and solve 
them numerically. The evaporation efficiency has a max- 
imum value where the advection-dominated structure of 
the corona (at large distances from the black hole) changes 
to a radiation-dominated structure (at small distances). 
This maximum can be related to the observed spectral 
transitions in black hole X-ray binaries. Our solutions of 
coronal structure are valid for both, accretion in galactic 
X-ray transients and in AGN. 

Key words: accretion disks - black hole physics - X-rays: 
stars - stars: individual: Cygnus X-1, Nova Muscae 1991 
- galaxies: nuclei 



1. Introduction 

Recently a wealth of X-ray observations has provided more 
detailed information about spectra of X-ray sources and 
short-time variability. This concerns soft X-ray transients 
as good examples of stellar accreting black holes, high- 
mass X-ray binaries with a black hole primary star and 
also Active Galactic Nuclei. The new observations put 
more stringent constraints on the modelling of accretion 
processes. In this connection the change from accretion in 
a standard cool disk to a hot coronal flow is an important 
feature. 

Different solutions of accretion on to compact objects 
had been investigated in the past. An account can be 
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found in the review by Narayan et al. (1998). The de- 
scription of an optically thick, geometrically thin disk, 
now often called "standard thin disk" , goes back to work 
of Shakura & Sunyaev (1973), Novikov & Thorne (1973), 
Lynden-Bell & Pringle (1974) and others. A much hotter 
self-consistent solution of optically thin gas was discov- 
ered by Shapiro, Lightman & Eardley (1976), but found 
to be thermally unstable by Piran (1978). A solution for 
even super-Eddington accretion rates in which the large 
optical depth of the inflowing gas prevents radiation to 
escape was suggested by Katz (1977), Begelman (1978), 
Abramowicz et al. (1988) and others. This is an optically 
thick advection-dominated accretion flow. The fourth so- 
lution concerns accretion at low sub-Eddington rate (pro- 
posed by Shapiro et al. 1976, Ichimaru 1977 and Rees et 
al. 1982). Detailed theoretical work on this solution by 
several authors has proven that this optically thin regime 
of solutions for low accretion rates is relevant for accretion 
in many cases (Narayan & Yi 1994, 1995a, b, Abramow- 
icz, Chen k Kato et al. 1995, Chen 1995; Chen et al. 1995, 
1997; Honma 1996; Narayan, Kato & Homna 1997a; Naka- 
mura et al. 1996, 1997). 

The model of an optically thin advection-dominated 
flow (ADAF) allows to explain the accretion in soft X-ray 
transients during quiescence. Assuming the old picture of 
a standard thin accretion disk reaching inward towards the 
black hole the amount of X-ray flux is clearly inconsistent 
with the accretion rate estimated from the outer optically 
thick disk (McCHntock et al. 1995). An interior ADAF 
allows to understand the observations but requires the 
change from a standard thin disk to the hot flow at a cer- 
tain transition radius. The black hole X-ray transients are 
a laboratory for accretion on to black holes. The physics of 
the accretion processes around supermassive black holes 
in AGN is similar, which makes the investigation of ADAF 
even more important and fruitful. 

The change from a geometrically thin disk to the 
hot flow is an essential feature in this overall picture of 
accretion. The "strong ADAF principle" , formulated by 
Narayan & Yi (1995b), suggests that, whenever the ac- 
creting gas has a choice between a standard thin disk and 
an ADAF, the ADAF conflguration is chosen. This rule 
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provides a transition radius r^mtr for a given accretion 
rate (see Fig. 6b in Esin et al. 1997). In a different ap- 
proach Honma (1996) considered the interface of a hot 
optically thin disk with a cool standard thin disk apply- 
ing radial diffusive transport of heat. 

Here we present a model of a corona above a thin disk 
in order to determine the inner edge of the thin disk, in- 
side of which only a hot optically thin flow exists. The co- 
existence of a cool disk and a corona had been investigated 
for disks in dwarf nova systems, where features otherwise 
not understood as for example the X-rays observed in qui- 
escence can then be explained (Meyer & Meyer-Hofmeister 
1994, Liu et al. 1995). The situation is the same in disks 
around black holes (Meyer 1999), both stellar and super- 
massive. 

The physical two-dimensional situation would demand 
the solution of a set of partial differential equations in ra- 
dial distance r and vertical height z. In particular a sonic 
transition would require treatment of a free boundary con- 
dition on an extended surface. In Sect. 2 we show how we 
derive a simplified set of four ordinary differential equa- 
tions which describe the physics of the corona above the 
thin disk. The interaction of corona and thin disk results 
in evaporation of matter from the thin disk . In Sect. 3 we 
discuss the general features of this evaporation process. 
The results of our computations are presented in Sect. 4. 
The main result of this model is the efficiency of evap- 
oration and with this the location where the thin disk 
has completely changed to a hot coronal flow (Sect. 5). 
In Sect. 6 we point out the invariance of our results for 
different black hole masses. Sects. 7 and 8 concern the va- 
lidity of the assumptions made in our computations and 
the importance of synchrotron and Compton processes. In 
Sect. 9 the comparison with observations follows, in Sect. 
10 a discussion of further effects. In Sect. 11 we compare 
with the model of Honma (1996). Sect. 12 contains our 
conclusions. 



2. The equations 

In this section we derive the equations that describe the 
corona above the inner rim of the accretion disk. For this 
we will average over an inner ring area and thereby reduce 
the set of partial differential equations to a system of or- 
dinary differential equations with respect to height. The 
solutions then become a function of the radial coordinate 
of the inner rim. 

We start with the five equations of viscous hydrody- 
namics, conservation of mass 



and the first law of thermodynamics 



| + V.(pv).0, 

the equations of motion 



(1) 



9v 



+ (v • V) V 



-VP + V-S- pV#, 



-P(V-v) + (S-V)-v + g, (3) 



supplemented by the equation of state for an ideal gas 

3? 



P 



-Tp, 



U = 



1 P 
7 - 1 p ' 



3? 



T = Vf 



(4) 
(5) 
(6) 



Here p, P and T are density, pressure and temperature, 
V is the velocity vector, S the viscous stress tensor, the 
gravitational potential, U the internal energy, 3? the gas 
constant, p the molecular weight. We take p = 0.62 for a 
fully ionized gas of cosmic abundances. Vg is the isothermal 
speed of sound, 7 the ratio of the specific heats, = |, for 
a mono-atomic gas. For S we take the form appropriate 
for a mono-atomic gas with Cartesian components 



O'ik = Pv 



dvj dvk_ _ 2 chn 
oxk oxi 3 oxi 



(7) 



(Sommerfeld 1949). For the (non- molecular) viscosity we 
use an a-parameterization (Shakura & Sunyaev 1973, 
Novikov & Thorne 1973) 



Mv = -^aP/f2, 



(8) 



fl rotational frequency, q is the rate of heating per unit 
volume. We write it in two parts, 

i = 5rad + 9cond 

where ^rad is the gain or loss of heat by radiation and ^cond 
that by heat conduction 



^cond — V • Fc 



(9) 



For (/rad we take here optically thin thermal radiation loss 
according to Fig. 1 of Raymond et al. (1976). 



?rad = -nenA(T) 



(10) 



where rie and n are electron and ion particle densities. For 
the heat flux Fc we take 

Fc = -KoT^/^VT; (11) 

(Spitzer 1962) for a fully ionized plasma and use the value 

^„_6 gcni 



neglecting a weak logarithmic dependence on electron 
(^-^ density and temperature. 
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Taking the scalar product of Eq.(2) with v, adding to 
Eq.(3), and using Eq.(l) we obtain the equation of con- 
servation of energy 



d_ 

dt 



■v[pU + p—+ p^ + P 



-(S-v)+Fc 



— Qrad 



(12) 



We introduce cylindrical coordinates, r, ^p, z with the z- 
axis perpendicular to the midplane of the disk and through 
the central accretor, and use the vector correspondence to 
the Cartesian advection term 



(v ■ V)v = Vy - V X (V X v), 



(13) 



where i)^ = (v • v) . The (Newtonian) potential is then 
GM 



$ = - 



(r2 + ^2)1/2 



(14) 



with G the gravitational constant and M the central mass. 

We consider from now on stationary, azimuthally sym- 
metric flows. 

In the following subsections we derive values for the 
rotational velocity v^p and the radial diffusive velocity Vr, 
and the three remaining ordinary differential equations 
for the conservation of mass, the vertical dynamic equilib- 
rium, and the energy flow in the corona above the inner 
zone of the accretion disk. 

For this we choose a region between inner disk radius 
ri and radius r2 such that 



7r(r| 



nr 



(15) 



where r is an area weighted mean radius between n and 
r2- We apply therefore, at each radius r, the integration 
bounds ri=0.72r and r2=1.23r. The coronal solution de- 
pends on the value of r as a measure of the distance from 
the accretor. 

We estimate the relative order of magnitude of the 
various terms in our equations and drop non-dominant 
terms. For these estimates we use relations that result 
from our solutions which are therefore self- consistent. 

All our solutions have 



GM 



GM 



<l/8, 



(16) 



see Fig. 3. They are subvirial and therefore like "slim" 
accretion disks (Abramowicz et al. 1988) with cylindrical 
symmetry in their (main) part in contrast to the more 
spherical symmetry of ADAF solutions. As a result 



Vr 



a 



(17) 



where we have replaced the azimiithal velocity by its Ke- 
pler value which is a good approximation in the main body 



of our solution, see the remark foUwing Eq. 25 below. A 
wind from such a corona reaches escape speed at large 
height leading to 



2r' 



(18) 



see Fig. 2. Where in the following we neglect radial deriva- 
tives of velocity components we use the order of magnitude 
relation 



or r 



2.1. The r-component of the equation of motion 



(19) 



This equation yields the rotational velocity. 

The radial component of the equation of motion be- 



comes 



p Vr- 



dVr 

dr 



dVr 

dz 



dP GMr 

~ 9^ ~^(r2 -hz2)3/2 

Id, . a 

+ - — {r(Jrr) + —C^rz 



r dr 



dz 



(20) 



In this equation the pressure term dP/ dr is small com- 
pared to the gravitational term in regions z < r/2 domi- 
nant with respect to dynamics and energetics of the solu- 
tions (see Fig. 2). Taking dP/dr =0{P/r) we estimate 



dP 
dr 



pGMi 



(r2 + ^2)3/2 



~ GM 

4^- 



(1 



)3/2 



r)^/^ < 0.18 



(21) 



This shows a relative unimportance of radial pressure 
support in these coronae (in contrast to virialized ADAF 
solutions, Narayan & Yi (1995b)). At a large height z = r 
this ratio has increased to 0.35. But we note that with the 
topology of Fig. 1 that at that height the radial pressure 
gradient must have significantly decreased because of the 
axial geometry. Therefore to first approximation we can 
neglect the pressure term. 

The viscous terms 



dr 



^V-v), arz 



,dVr 

dz 



dvz 
dr 



(22) 



contribute derivatives of r- and ^;-components of the veloc- 
ity. Their contributions, evaluated with Eqs. (16) to (19), 
are small compared to 



dP 
dr 



0{P/r) 



(23) 



by at least a factor Vs/v^. They as well as the first two 
terms on the left hand side of Eq. (20) are negligible. 
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This leaves the balance between centrifugal and radial 

gravitational force 



n = 



GM 



GM 

r3 ^^'^ r^' 



:L -1-3/4 



(24) 



(25) 



In the main body of the coronal solution the Kepler 
values (i.e. at z^Q) are a good approximation to v^p and 
i7. In our calculation we take the Kepler values throughout 
as a computational convenience. It overestimates the fric- 
tional work and underestimates the release of rotational 
energy at large height. These two effects partially com- 
pensate each other. At z < r/2 the sum of these effects 
deviates from the accurate value by less than 7%. At z = r 
it is 20%. But at this height the total further contribution 
of both energy terms already is negligibly small (due to the 
exponential decrease of pressure and density with height, 
see Fig. 2) and isothermality of the solution completely 
dominates. This simplification has thus only a small effect 
on the accuracy of the solution. 

2.2. The if-component of the equation of motion 

The conservation of the z-component of angular momen- 
tum yields the r-component of the velocity. 

Multiplication of the equation of motion by 
rUip{ur,u^,Uz unit vectors in r, (p and z direction) 
yields 



pv ■ V{rv^) = ruy • (V ■ 5) 

= V • (S • ru^) - (S • V)ru^ (26) 

Using (1) (for f =0) 

V ■ pv = 0, 
and 

(the right side understood as the dyadic product) and 

this results in the equation of conservation of 2-angular 
momentum 

V • {pwrv^ - rvi^ • S) = 0. (27) 

We integrate over a cylinder with radius ri and obtain 



27rri^ / pVrVtpdz = 2-Kr\ 
Jo 



/>oo 

/ Uripdz. 

Jo 



(28) 



Here it is assumed that no significant angular momen- 
tum is accreted at the center of the compact object or 
lost through the boundary z = 00. For accretion and jet 



angular momentum loss from the most interior parts of 
the accretion region, such losses would be small by order 
of ^JR|r\ where R is the radius of the inner accretion 
region, e.g. R « 3rs for the Schwarzschild radius rg of a 
central black hole (rs = 2GM/c^). 

We replace the global balance of angular momentum 
flow of Eq. (28) by its local equivalent and obtain the 
radial diffusive velocity 



PVip 



-a- 



(29) 



with ffrip = pv{dvif,/dr ~ v^p/r) = —aP. This is the stan- 
dard bulk radial velocity formula for accretion disks. It 
neglects the modification by ^-dependent secondary flows 
superimposed on the main drift of material (see Kley & 
Lin 1992). Such secondary flows may slightly modify the 
local vertical structure but will hardly affect the general 
structure of the solutions. 



2.3. The z-component of the equation of motion 

The vertical component of the equation of motion Eq. (2) 
yields the vertical dynamic equilibrium. 
The vertical component becomes 



dv^ dv^.__dP_ GMz 



Id, 
r or 



d 



(30) 



Using approximate isothermal pressure layering (see 
Fig. 2) 



(31) 



with hK.\r (see Fig. 3) and the estimates of Eqs. (16) to 
(19) we estimate 



(32) 







dP 


or 


1 


dz 



< 0.02a 

to be negligible. Further, with 

dvr dv. 



(Jrz = Py \ -5 1- 

oz or 



^- = /^vl2^-|(V-v))(33) 



the largest contribution of all viscous terms in Eq. (30) 
appears in ddzz/dz) from the z-derivative of p^ (that is 
of P, see Eq. (8) times the derivative dvz/dz. This term 
is of order <0.16a RiO.05 and negligible. In the remaining 
dominant terms only z-derivatives appear. Averaging over 
the inner zone we obtain 



PVz 



dvz 
dz 



__dP_^ GMz 



dz (r2 + 02)3/2- 



(34) 
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2.4- Mass conservation 



2.5. Energy conservation 



The equation of mass conservation we integrate over the 
radial zone between ri and r^ , of area 7rr^ and obtain 



We first derive the frictional fliix terms. With Eqs. (22) 
and (13) we obtain 



r2 Q 

2'jTr—(pVz)dr 
oz 



{2'KrpVr)l — {2'KrpVr)2- 



(35) 



When, as we 
strongly (e.g. p ~ 



assume, rpv^. drops off with radius 
r~^ in the scaled solutions, Liu et al. 
1995) the difference on the right hand side is dominated 
by the first term, and its value will be slightly smaller 
than this term. We approximately take this into account 
by replacing the term at ri by its value at the somewhat 
larger mean radius r. We replace the left hand side by the 
area times a mean value. This leads to the approximation 



A. 

dz 



{pVz) 



-pVr 



(36) 



(v ■ S) = Mv 
and its r-component 



Vw^ - V X (V X v) - -v(V • v) 



Ur • (v • S) = Py 



dv 



dr dr r 

2 



dvr dvz 
dz dr 



(40) 



(41) 



Here all terms are negligible except those containing 
v^p, by use of Eqs. (16) to (19). Making use of Eqs. (8) 
and (19) we obtain 



where is the radial flow velocity (Eq. 29). At large 
height, however, we have to take into account the flaring 
of the channel as the ascending gas changes from a cylin- 
drical rise to spherical expansion (see Fig. 1). How can we 
account for the additional attenuation of the vertical flux 
density by this change of geometry? 

We proceed in the following way. Approximately the 
cross section of the vertical flow channel can be taken as 
proportional to 1 + z^/r^. Without sidewise diffusive loss, 
mass conservation would then yield 



Ur • (v ■ S) = —aPv^ 



(42) 



Likewise, only keeping non- negligible v^p— terms, the z- 
component is 



• (V • S) = /Zv 



dz 



Vr 



dvr dv. 



dz 



dr 

2 dv^ 



(43) 



d_ 

dz 



1 + 12 M« 



= 



The radial component of the divergence can be written, 
(37) using (19) as 



if we take the center of the expanding flux tube in the 
vertical direction. Per unit cross section the resulting con- 
servation equation is then 



l + 4dz 



pVz 



d 

dz 
= 



pVz 



2z 

H + Z'' 



(38) 



The gas in a cylindrical column thus experiences an 
additional geometry-related radial divergence 



2z 



(r^ + z2) 



PVz 



Adding the two radial flow components we obtain the final 
form 



d , , 2 2z 



(39) 



We note that the diffusive part in this equation be- 
comes unimportant at large height where the diffusive de- 
scription for the flaring flux tube becomes problematic. 



-aPn\ - + = ^aPf}. (44) 

, 2 alnr / 2 



For 91nP/i91nr we have here taken its estimated value 
-2 at the characteristic height z = h = r/2 (see Fig. 2), 
which we obtain from Eq. (31) with the scaling Pq ~ ''"^ 
that holds when we compare solutions with different dis- 
tances r (see Liu et al. 1995). The resulting term repre- 
sents the dominant heating mechanism of the corona by 
release of gravitational energy. 

We now integrate over the area between ri and r2 as 
in the derivation of conservation of mass 



zirrdri — 
[dz 



Vz[P^ 



7-1 



P + p<P 



(45) 



= i 27rr 



Vr[PY + :^P + P^]+{Fo)r 
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On the right hand side of this equation we again re- 
place the difference of the terms taken at ri and r2 by 
the value taken at the intermediate value r, because the 
value at r2 will only be a fraction of the value at ri due 
to the radial drop off of all the quantities. If we take Vr as 
the diffusive velocity (Eq. (29) we must still include the 
divergence of the vertical flow due to the flaring geometry. 

We use the same procedure as in the preceding Sect. 2.4 = [Fc)z = —koT^^"^ — 
to derive the corresponding term and add it on the right 
hand side. Again, on the left hand side of the equation we 
write the integral as area times a mean value and obtain 
the energy equation in its final form 



This follows when one uses the scaling of P, T and f2 
derived by Liu et al. (1995) (see also values Pq in Table 
1). In the "conduction-dominated" corona r < r (Merit) 
the temperature becomes constant and the ratio becomes 
zero. The contribution of (i^c)r can therefore be neglected. 
The remaining ^—component of Eq. (11) we write as 



dT 



dz 



(49) 



'. 6. The system of ordinary differential equations with 
boundary conditions 



d_ 

dz 



\aPO + Qrad 



2 

+ - 

r 



P + p$\+iF,)r 



2z 



r^ + z^ 



(46) 



Had we chosen to integrate the radial divergence 
[rUr • (v • S)] underneath the integral of Eq. (45) in- 
stead of evaluating this term as in Eq. (44) the last equa- 
tion would have the factor 2 instead of |. This difference 
arises from the different assumptions about radial varia- 
tion of the quantities, implicit in the two approaches. It 
shows a natural limit of numerical accuracy of any such 
approach. 

The lateral inflow of heat by (Fc)r adds to the fric- 

tional release of energy in the coronal column. We show 
that the relative size of the contribution is small and we 
will take (Fc)r=0 in our computations. We assume that 
the temperature at low z/r varies like the virial tempera- 
ture, T ~ 1/r, and that at z > r the radial temperature 
gradient becomes small as the topology changes and the 
corona fills the full cylindrical region, with ^ = on the 
axis of symmetry. With these assumptions we take 



dT 
dr 



T 



(47) 



The ratio of total radial conductive heat input {Qc)r to 
the total frictional energy input Qf is then evaluated as 



{Qc)r 

Qf 



ld_ 
r dr 



[r{Fe)r] dz 



-aPQdz = 0.22, 



(48) 



It is independent of the value of r as long as we are in 
the "advection-dominated" ranger > r(Mc.rit) (see Fig. 3). 



Prom the five partial differential equations (1) to (3) we 
have now derived three ordinary differential equations 
(34), (39) and (46) (where we take (Fc)r = 0) which de- 
scribe vertical dynamical equilibrium and conservation of 
mass and energy, plus the two expressions for and Vr 
(Eqs. (24) and (29)). Our fourth ordinary differential equa- 
tion is the dependence of the heat flux on the temperature 
gradient, Eq. (49). As the 4 dependent variables we have 
temperature T, pressure P, mass flow m = pvz and heat 
flux Fc- The 4 boundary conditions are: 

(1) No pressure at infinity (i.e. no artificial confine- 
ment). This requires sound transition at some height 
z = Zi (free boundary problem) 



Vg a.t z = z\. 



(50) 



(2) No influx of heat from infinity. We neglect a small 
outward conductive heat fltix induced by the temperature 
gradient in the expanding wind and require 

Pc = at z = zx. (51) 

(3) Chromosphere temperature at the bottom z = z^ 



T = Tcff at z = zq. 

(4) No heat inflow at the bottom 

Fc = at Z — Zq. 



(52) 



(53) 



This approximates the condition that the chromospheric 
temperature is kept up by other processes and cannot con- 
duct any sizeable thermal heat flow. 

The consequence of conditions (3) and (4) is that the 
conductive flux that comes down at a certain level above 
the boundary must be radiated away. Since the advective 
cooling by the rising mass flow is small (in our solutions 
at T^-^K only ^) this reduces to the balance between 
thermal conductive flux and radiation and establishes the 
temperature profile as the one studied by Shmeleva & Sy- 
rovatskii (1973) for the solar corona. In these lowest layers 
the pressure is nearly constant and the temperature rises 
steeply. Liu et al. (1995) recalculated the temperature pro- 
file using the more recent cooling function of Raymond et 
al. (1976) and we used the improved version. This allows 
to replace the conditions (3) and (4) by 

r = 10S-^K and Fc = -2.7310Sp at z = zq (54) 
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The relation between P and Fc is practically independent 
of the value of the chromospheric temperature as long as 
it is small compared to Tq = 

These equations and boundary conditions describe a 
model of the corona-disk interaction for the innermost re- 
gion of the thin accretion disk where evaporation is most 
efficient. 



in the corona as a consequence of conservation of angular 
momentum. (One might compare this with the flow in a 
"free" thin disk without the tidal forces acting in a binary. 
In such a disk matter flows inward in the inner region and 
outward in the outer region, with conservation of the total 
mass and angular momentum (Pringle 1981)). 



3. General features of evaporation 

The equations derived in the last section describe the 
physics of a hot corona in equilibrium with a thin cool 
disk underneath. Aspects of this model had already been 
investigated for accretion on to a white dwarf (Meyer & 
Meyer-Hofmeister 1994, Liu et al. 1995, for application to 
WZ Sge stars and X-ray transients see also Mineshige et al. 
1998). The picture that arises from the processes involved 
is the following. Heat released by friction in the corona 
flows down into cooler and denser transition layers. There 
it is radiated away if the density is sufficiently high. If the 
density is too low, cool matter is heated up and evapo- 
rated into the corona until an equilibrium density is estab- 
lished. Mass drained from the corona by an inward drift is 
steadily replaced by mass evaporating from the thin disk 
as the system reaches a stationary state. This evaporation 
rate increases steeply with decreasing distance from the 
central compact star. When the evaporation rate exceeds 
the mass flow rate in the cool disk, a hole forms up to the 
distance where both rates are equal. Inside this transition 
radius rtr only a hot coronal flow exists. 

The mass flow in the corona is similar to that in the 
thin disk. The matter has still its angular momentum, 
and differential (Kepler-like) rotation causes friction. The 
main part of the coronal gas flows towards the center di- 
rectly. For conservation of angular momentum a small part 
flows outward in the coronal layer and condensates in the 
outer cool disk. Since the temperature is very high the 
corona is geometrically much thicker than the disk un- 
derneath. Therefore processes as sidewise energy trans- 
port and downward heat conduction are important for the 
equilibrium between corona and thin disk, together with 
radiation and wind loss. 

We model the equilibrium between corona and thin 
disk in a simplified way using a one-zone model as de- 
scribed in Sect. 2. This is possible since the evaporation 
process is concentrated near the inner edge of the thin disk 
and we can choose a representative radial region from ri 
to r2 such that the main interaction between cool disk 
and corona occurs there and evaporation further outward 
is not important. In Fig. 1 we show the flow pattern as 
derived from analysis and simplified modelling for the case 
that a hole in the thin disk exists. There are three regimes. 
(1) Near the inner edge of the thin disk the gas flows to- 
wards the black hole. (2) At larger r wind loss is impor- 
tant taking away about 20% of the total matter inflow. 
(3) At even larger distances some matter flows outward 




outer thin accretion I one zone_^ 
disl< model region 



Fig. 1. Coronal mass flow pattern (schematic). Convective 
circulations will be superimposed. 



4. Results from computations 

4-.1. Computer code technique 

We solve the differential equations (34), (39), (46), and 
(49) using the Runge-Kutta method integrating from the 
lower boundary at an assumed height zq ^ r above the 
midplane to the upper boundary at zi . The value Zi is de- 
termined by our computations as the height where sound 
transition occurs (free boundary problem) . The boundary 
conditions are given in Eqs. (50), (51) and (54). We re- 
quire that the energy flux at the upper boundary is zero 
but are content in our computations if it is very small, 
< 0.005 •|(F,)^ax|. 

We have to vary Pq and toq, trial values for P and 
m (m = pvz) ai z — zq, to find the coronal structure 
which fulfills the upper boundary conditions. The condi- 
tions, sound transition and zero heat flux, are met only for 
a unique pair of values for pressure and mass flow at the 
bottom. For other values sound transition is not reached 
and/or the heat flux is not neghgible (this is the same situ- 
ation as in coronae of cataclysmic variable disks, compare 
Fig. 2 of Meyer & Meyer-Hofmeister 1994). Numerically, 
the smaller the distance r from the compact object the 
more accurately the pressure value has to be determined 
to obtain a consistent solution. 

For our computations we put the lower boundary at 
about Zq — 2 10® cm above the midplane, but for r < 3 10* 
cm at correspondingly smaller values zq. In most cases the 
true height of the thin disk is less than the chosen value zq. 
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We carried out test computations to check the influence 
of assuming different values zq on the coronal structure. 
As expected the changes were found to be very small. 

Pq and (-F'c)o are nearly independent of the value of the 
chromospheric temperature as long as it is small compared 
to To = 

Since our results also concern accretion in AGN we 
have confirmed the validity of our assumption for the lo- 
cation of the lower boundary for this case also. 

4-2. Coronal structure 

We took 6Mq as our standard parameter value to model 
the evaporation in accretion disks around stellar black 
holes. We compute the structure of the corona for dis- 
tances of the outer cool disk from the black hole rang- 
ing from large values to small values of a few hundred 
Schwarzschild radii. For different distances r the struc- 
ture in z direction is similar if considered as a function of 
z/r. As will be discussed in Sect. 6 the coronal structure 
scales with Schwarzschild radius and Eddington accretion 
rate. 
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Fig. 2. Coronal structure at distance r = 10^'^cm from 
the 6Mq black hole. Temperature T is given in units of 
0.2-virial temperature, pressure P and vertical mass flow 
p ■ Vz (vz vertical velocity) are scaled to their values at the 
lower boundary. -Fc is the downward conductive heat flux, 
measured in units of the frictional energy generation rate 
in a column of unit cross section. Vg is the sound velocity. 

In Fig. 2 we show the structure at r = 10^'^cm close to 
where the evaporation efficiency is highest for a 6Mq black 
hole. Pressure and density decrease significantly within a 
narrow range above the cool disk. The downward conduc- 
tive heat flux -F^ therefore also has its maximum at a low 



height z, then rapidly decreases with height and is neg- 
liglibly small where sound transition Vz/Vg = 1 is reached, 
our boundary conditions. This happens at 2.3 lO^cm. The 
temperature is practically constant in the higher layers 
of the corona, a consequence of the high thermal conduc- 
tivity. The computed coronal structure at distance r cor- 
responds to the situation that the cool disk underneath 
extends from this distance r (inner edge) outwards. The 
derivation in Sect. 2 is based on the dominant role of the 
innermost corona and our model cannot be used to de- 
scribe the coronal structure above the cool thin disk far- 
ther out which will be affected by the dominant process 
farther in. 

These results can be compared to those of our earlier 
investigation of evaporation in cataclysmic variable disks. 
The coronal structure was considered there at distances 
10* '' to 10^ '' cm from the IMq white dwarf adequate for 
a disk around a white dwarf. Due to the scaling with the 
mass of the compact star the structure above the thin disk 
around a IMq white dwarf at distances logr 9.0 and 9.7 
(shown in Liu et al. 1995, Fig. la and lb) corresponds to 
the structure at distances logr/rs = 3.53 and 4.23 from a 
6Mq black hole. 
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Fig. 3. Values of quantities at the inner edge r=rtr of 
the thin disk for various distances r (in cm). Rate of in- 
ward mass flow Af (in g/s) in the corona {— evaporation 
minus wind loss), maximum temperature in the corona 
and h/r {h pressure scaleheight), solid lines. Dashed line: 
virial temperature. Values of M and r are for a central 
mass of 6Mq, scaling with MEdd and rg is indicated. 

With distance from the black hole the coronal struc- 
ture changes. The main result concerns the evaporation 
rate, the shape of the function M. For small values of r 
the evaporation rate M increases with r, then reaches a 
maximum value and decreases hirthcr outward. In Fig. 3 
we show the evaporation rate, the temperature reached at 
sound transition and the ratio of the pressure scaleheight h 
to r (maximal temperature and pressure scaleheight from 
the computed vertical structure of the corona). In Table 
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Note: rho and Pq vertical mass flow rate and pressure at the lower boundary zq, T temperature reached at the upper boundary 
(sound transition), M rate of inward mass flow in the corona (=evaporation rate); A fraction of mass carried away by the wind; 
quantities r/rs and M/MEdd scaled to Schwarzschild radius and Eddington accretion rate, r, mo,Po in cgs units. 



1 we summarize values of characteristic quantities. Note 
the large variation in the values mp and Pq: the values 
of P and rh at the lower boundary. The evaporation rate 
is given in {M(7)/y) and scaled to the Eddington accre- 
tion rate (MEdd — iEdd/O.lc^), the radius in units of 
Schwarzschild radius r/rg is added. For different black 
hole masses, the evaporation rate can be scaled accord- 
ingly (Sect. 6). 

To include the evaporation process in the evolution of 
the accretion disk, e.g. for black hole X-ray transients, the 
M-r relation for values r larger than the values in Fig. 3 
might be needed. For such distances much larger than the 
hump location (Fig. 3, logr/rs > 3) the function AI can 
be approximated as found in the earlier investigation by 
Liu et al. (1995) for accretion in cataclysmic variables. 

In Table 1 we also give the fraction of mass carried 
away by the wind. Note that this fraction is of the or- 
der 20% at large distances r and decreases to very small 
values inside the maximum of the function M (Fig. 3). 
The location of sound transition is about at the height 
3r for a distance r — 10^'^cm, the hump location, and at 
increasingly larger heights further inward. But their up- 
permost layers are unimportant for the energy balance in 
the corona because of their very low density. 

Disks around black holes can reach inward even to 
the last stable orbit, and a new feature appears, a maxi- 
mum efficiency of evaporation. In dwarf nova systems the 
inner disk edge disk lies much further out measured in 
Schwarzschild radii. The maximum of the evaporation ef- 
ficiency has important consequences for the inward ex- 
tension of the thin disk (see Sect. 5). This maximum is 
caused by the change in the physical process that removes 
the heat released by friction. For larger radii the frictional 
coronal heating is balanced by inward advection and wind 
loss. This fixes the coronal temperature T at about 1/8 of 
the virial temperature Tv {^kT^/ ji = GM/r). Downward 
heat conduction and subsequent radiation in the denser 
lower region play a minor role for the energy loss though 



they always establish the equilibrium density in the corona 
above the disk. The ratio of the pressure scaleheight h to 
r is almost constant for larger r and becomes small closer 
to the compact object. 



4-. 3. The energy balance in the corona 
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Fig. 4. Contributions to the local energy balance in the 
corona above the thin disk around a 6Mq black hole at 
logr=8.8. Solid lines: gain by friction and vertical advec- 
tive fiow, dashed lines: loss by radiation, vertical heat flow 
and sidewise advective outflow (flow out of the unit vol- 
ume) 



10 



F. Meyer, B.F. Liu, E. Meyer-Hofmeister: Evaporation from accretion disk to coronal fiow 



To illustrate the local energy balance we show in Fig. 4 
all contributions separately as functions of height z. There 
are two terms of energy gain, by friction and by verti- 
cal advective flow. Friction decreases with height with the 
pressure (compare the coronal structure in Fig. 2). The 
gain by vertical advection results from the fact that more 
mass enters from the layer below than leaves into the 
layer above. This is due to the net sidewise outflow out 
of the unit volume. This overcompensates the increase of 
enthalpy with height. Concerning the losses, the radiation 
loss is the essential part in lower layers. It is proportional 
to the square of density and decreases steeply towards 
higher z. The loss by vertical heat flow describes the ef- 
fect that the heat, released by friction and advection in 
higher layers is conducted down to layers where it is radi- 
ated away. 

When the temperature has dropped to at very 

low height the downward heat flow is then 10 times larger 
than the upward advective flow. The balance is then be- 
tween conduction and radiation. This leads to the temper- 
ature proflle calculated by Shmeleva & Syrovatskii (1973) 
and Liu et al. (1995) and fixes the relation between and 
P at this temperature used as boundary condition for our 
calculations. 

For logr=8.8 (Fig. 4) the height-integrated losses by 
advective outflow and by radiation are about equal. This 
is a particular case. With distance r the relative size of 
these two losses varies in a characteristic way. To show 
this we have integrated the two losses over the height of 
the (flaring) column and normalized by the total release of 
energy by friction. Fig. 5 shows the result. The advective 
loss term includes the ernergy carried advectively away 
with the wind. This latter contribution is also shown sep- 
arately. 

The maximal coronal mass flow rate is reached at that 
radius where the character of the solution changes from 
an advection-dominated flow at large radii to a radiation- 
dominated flow at small radii (see Fig. 3). (This is an 
intrinsic feature of the solutions as a dimensional analysis 
of the global energy balance equation can show). 

This behaviour contrasts with that of optically thin 
accretion flows of constant mass flow rate M. These 
are advection-dominated (ADAF) at small radii but be- 
come radiation-dominated (Shapiro, Lightman and Eard- 
ley 1976) at large radii (Narayan & Yi 1995b). The coronal 
evaporation, on the other hand, allows for mass exchange 
between disk and corona and results in a M varying with 
r which results in an inversion of that behaviour. 

This mass exchange also provides thermal stability to 
the coronal evaporation solution. The system responds to 
a slight increase in temperature by an increased thermal 
flux that leads to an increase of density due to increased 
evaporation and thereby to a decrease of temperature by 
increased radiative cooling. 
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Fig. 5. Distribution of the total energy released by fric- 
tion. Solid lines: fraction going into radiation and fraction 
removed by advection (sidewise and via the wind); dashed 
line: wind part only. (6M0 black hole) 



5. The inner edge of the thin disk and spectral 
transitions 

If the mass flow rate in the disk is low, the evaporation 
rate in the inner disk region may exceed the mass flow 
rate. The inner edge of the cool disk then is located where 
both rates are equal, that is where all matter brought in 
by the thin disk is evaporated. Inside this location only a 
coronal flow exists. This coronal flow provides the supply 
for the ADAF towards the compact object. The higher the 
mass flow rate in the cool disk the further inward is the 
transition to a coronal flow/ADAF. This is the situation 
in transient black hole X-ray binaries during quiescence. 
If the mass flow rate exceeds the maximal evaporation 
rate, M > Merit (the peak coronal mass flow rate), the in- 
ner cool disk cannot be completely evaporated anymore, 
instead it reaches inward to the last stable orbit. In high- 
mass X-ray binaries fluctuations of the mass flow rate have 
the same effect. In Fig. 6 we show how the inward extent of 
the standard thin disk depends on the mass flow rate in the 
thin disk. As indicated there the corresponding spectra are 
hard, if the hot coronal flow/ADAF dominates, or, soft, if 
the standard thin disk radiation dominates near the cen- 
ter. This picture from our modelling is the same as the de- 
scription of the accretion flow for different spectral states 
presented in the investigation of advection-dominated ac- 
cretion by Esin et al. (1997, Fig. 1, except the very high 
state). 

As illustrated in Fig. 5 the transition between dom- 
inant advective losses further out and dominant radia- 
tive losses further in causes the change in the slope of 
the evaporation rate at M= Merit- Except for the differ- 
ence between the sub-virial temperature of the corona and 
the closer-to-virial temperature of an ADAF of the same 
mass flow rate, a similar critical radius is then predicted 
by the "strong ADAF proposal" (Narayan & Yi 1995b). 
In general, however, the strong ADAF proposal results in 
an ADAF region larger than that which the evaporation 
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model yields. A detailed discussion of these spectral tran- 
sitions was given separately (Meyer et al. 2000). 
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Fig. 6. Evaporation rate M{r) as in Fig. 3. Inward ex- 
tension of the standard thin disk for 3 different mass flow 
rates M in the thin disk (schematic). Note that the stan- 
dard thin disk reaches inward towards the black hole if 
M > -Merit- Shown is also the type of spectrum, soft or 
hard, related to M 



The main results from our model are the following. 
The transition from a thin disk to an ADAF should oc- 
cur at about the same luminosity for different systems if 
their black hole masses are not too different. The inner 
edge location should then either be at a distance of some 
Schwarzschild radii (the distance for which M = Merit) or 
the thin disk should reach inward to the last stable orbit. 
Several high-mass binaries and black hole transients show 
such a systematic behaviour (Tanaka & Shibazaki 1996, 
Tanaka 1999). But in some cases the transition seems to 
be farther inward. In Sect. 9 we compare our results with 
observations. 

6. Solutions for different central masses - 
stellar and supermassive black holes 

Our equations and boundary conditions are invariant 
against scale transformations that leave temperature and 
velocities constant. An invariant Kepler velocity implies 
that M/r is invariant as is Fi.r,Pr,pr,mr and M/r. One 
thus obtains solutions for another central mass M2 from 
those for Mi by multiplying r, z and M values by M2/M1, 
and by dividing Fc,P,p and to values by the same factor 
while keeping velocities and temperatures the same. Since 

AQttGM 2GM 
MEdd = , rs = — 5— (55) 

KC & 

are both proportional to mass, the relation between M and 
r is invariant if plotted in units of Eddington accretion rate 



and Schwarzschild radius. In these units Fig. 3 is universal 
and applies to both stellar size and supermassive black 
holes. 

The properties of the thin disk underneath do not scale 
similarly. This has no effect on the coronal structure, but 
shows up in the spectra belonging to a thin disk -|- corona 
for a given mass accretion rate. 

7. Validity of assumptions 

1.1. Heat conduction 

In our modelling classical heat conduction is assumed. It is 
only valid if the conductive flux remains sufficiently small 
against the maximum transport by free streaming elec- 
trons. The relevant ratio is estimated here as 0.16 at the 
radius of maximal mass flow, constant for smaller radii, 
and decreasing as l/r~^/^ for large radii. All estimates 
are independent of the central mass. Smaller a-values in- 
crease the validity of the approximations. 

1.2. Equilibrium between electrons and ions 

Our model assumes a one-temperature plasma. Tempera- 
ture equilibrium between electrons and ions requires that 
the relevant collision timescale between both species re- 
main small compared to the heating timescale. Their ra- 
tio is estimated as 0.27 where the coronal mass flow rate 
reaches its maximum value, and scales as r~^/^. Thus at 
smaller radii this approximation brakes down and elec- 
trons and ions will develop different temperatures. 

8. Synchrotron radiation, Compton effect 

We have performed estimates on the size of these neglected 
effects. Synchrotron losses are proportional to the number 
density of electrons, their kinetic energy, and the square of 
the magnetic field strength. It is often assumed that the 
magnetic field is dynamo produced and that its energy 
density is a fraction ^ of the gas pressure. The ratio of 
synchrotron losses to thermal plasma radiation then be- 
comes a pure function of temperature rising with T3/2_ 
The main cooling in our solutions occurs in lower lay- 
ers at about 1/2 of the final coronal temperature. Un- 
der such circumstances, synchrotron radiation remains less 
than bremsstrahlung for /3 > 4. Values for /3 in the litera- 
ture include /3 > 5 from energy considerations for ADAFs 
(Quataert & Narayan 1999b) and /3 < 10 from direct 
numerical MHD calculations (Matsumoto 1999, Hawley 
2000). 

Compton cooling by photons from the underlying disk 
depends on the mass flow rate in the disk. For a rate 
equal to the maximal coronal mass flow rate, this Comp- 
ton cooling approaches about 1/3 of the frictional heating. 
Compton cooling (and heating) of the corona by radia- 
tion from the accretion center depends on the state. In 
the hard state, taking the example of Cyg X-1 (Gilfanov 
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et al. 2000), the spectrum yields a Conipton equilibrium 
temperature of 10*'^ K. From the flux one estimates a 
Compton energy exchange rate comparable with the fric- 
tional heating rate in the corona. This indicates that the 
true coronal temperature at the peak of the coronal mass 
flow rate of Fig. 3 could be somewhat smaller. A similar 
result is obtained for Compton cooling by the central flux 
during the soft state. The keV-photons efliciently cool the 
corona heated by friction when the temperature obtained 
with the pure bremsstrahlung assumption is larger than 
an estimated lO^ '^K. 

All these cooling mechanism tend to reduce the peak 
temperature and mass flow rate in Fig. 3 and shift the 
minimal transition radius somewhat outwards. Our esti- 
mates are only rough and a more quantitative evaluation 
of these eff'ects has to be done later. 

9. Comparison with observations 

Our model predicts a relation between the evaporation 
rate and the inner disk location (function M in Fig. 3 and 
6). This fact can be compared with observations in two 
ways: [1] for low mass flow rates in the thin disk a hole 
should form; the location of the change from accretion 
via the thin disk to a coronal flow should be related to 
the mass flow rate as predicted; [2] the spectral changes 
hard/soft and soft /hard should occur for a predicted mass 
flow rate and the location of the inner edge at this time 
should agree with that from the model. To perform test 
[1] we need to know two quantities, the mass flow rate in 
the thin disk and the inner disk edge location. Fitting of 
the spectrum based on the ADAF model (for a review see 
Narayan et al. 1998) gives values of the mass accretion rate 
on to the black hole and therefore also the evaporation 
rate (= mass flow rate in the disk at the inner edge). 
For the location of the inner edge constraints come from 
the Hq, emission line. The maximum velocity indicates the 
distance of the inner disk edge from the black hole. For 
test [2] the luminosity for which the spectral transition 
happens has to be considered, as well as signatures of the 
inner disk edge. 

9.1. Application to X-ray novae 

For the X-ray transients A0620-00 and V404 Cyg Narayan 

et al. (1996,1997a) derived mass accretion rates based on 
spectral fitting of the ADAF model. These rates were 
used to perform test [1] . For A0620-00 the resulting inner 
edge location agrees reasonably with the value from the 
Ra emission line. For V404 Cyg a discrepancy was found 
(Meyer 1999), but this disappears with the recent results 
for V404 Cyg (Quataert & Narayan 1999b) and we find 
agreement also for this system. The comparison for GRO 
J1655-40 using the rate from Hameury et al. (1997) and 
outer-disk-stability arguments is satisfactory. For details 
of these comparisons see Liu et al. (1999). 



The agreement of results derived from the disk evapo- 
ration on one side and the ADAF spectral fitting, includ- 
ing different physics for the description of the innermost 
region, on the other side, supports both models. 

Further agreement (related to test [1]) between ob- 
servation and theoretical modelling of evaporation comes 
from the computed disk evolution of A0620-00 (Meycr- 
Hofmeister & Meyer 1999). The long lasting outburst 
makes plausible that only little matter is left over in the 
disk. If this is the case the evolution during quiescence 
can be considered independent of the detailed outburst 
behaviour. We took into account the evolution of the thin 
disk together with the corona. The evaporation determines 
the inner edge of the thin disk at each time. The mass flow 
rate there is important for the disk evolution and therefore 
for the outburst cycle (without careful treatment of the 
inner edge boundary the evolution cannot be determined 
adequately). The evolution based on the same frictional 
parameter as taken for dwarf nova modelling (acooi=0.05). 
leads to the onset of an outburst after the right time and 
the amount of matter stored in the disk in agreement with 
the energy estimated from the lightcurve. 

Concerning test [2] spectral transitions observed for a 
few transients can be considered. The best observed source 
is Nova Muscac 1991 (Cui ct al. 1997). The transition oc- 
curs around Lx « 10'^'' erg/s (Ebisawa et al. 1994). These 
spectral transitions were modelled based on an ADAF by 
Esin et al. (1997). Our value for the critical mass accre- 
tion rate for a 6 M© black hole, 10^^'^g/s, corresponds to 
a standard accretion disk luminosity of about 10^''-^ erg/s, 
in agreement with observations. 

9.2. Application to high-mass X-ray binaries 

The accretion rates observed for high-mass X-ray bina- 
ries are higher than those for transients, fluctuations are 
common. Therefore hard/soft and soft/hard spectral tran- 
sitions occur much more often than for transients where 
this phenomenon is only related to outburst onset or de- 
cline. In our investigation of spectral transitions (Meyer et 
al. 2000) we discuss the behaviour of the three persistent 
(high-mass) black hole X-ray sources LMC X-1, LMC X-3 
and Cyg X-1. The transition occurs around Lx ~ 10"^^ 
erg/s as also for transients (Tanaka 1999). For the spec- 
tral transitions of Cyg X-1 again the ADAF model was 
successfully used to describe the observations (Esin et al. 
1998). 

Our model predicts the inner edge of the thin disk at 

some 100 Schwarzschild radii for the accretion rate related 
to the hump in Fig. 6, that is at the moment of spectral 
transition. We have estimated the timescale on which a 
standard accretion disk at accretion rates near the crit- 
ical value Merit and at the related radius can evaporate 
or form newly and find a value of the order of days. This 
agrees with the timescale for the observed spectral tran- 
sition of a few days for Cyg X-1 (Zhang et al. 1997). Fre- 
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quency resolved spectroscopy of Cyg X-1 (Rcvnivtsev et 
al. 1999) also leads to an inner edge of the thin disk at 
150 Schwarzschild radii. But the existence of a reflecting 
component indicates a disk further inward (Gilfanov ct 
al. 1998, Zycki et al. 1999). For a sample of X-ray tran- 
sients and high-mass X-ray binaries Di Matteo & Psaltis 
(1999) found, if they relate the rapid variability properties 
to those of neutron stars, that the disk reaches down to 
about 10 Schwarzschild radii. This would be in contrast to 
the prediction from the ADAF based modelling and our 
evaporation model. Further work is needed to understand 
these signatures of a thin disk so far inward during the 
hard state. 

9.3. Application to AGN 

Liu et al. (1999) used the results derived in the present pa- 
per to predict the location of the inner edge of the thin disk 
around M87 (test [1]). Using the accretion rate derived by 
(Reynolds et al. 1996) and an inner edge of the thin disk 
from emission lines broadening (Harms et al. 1994) they 
found good agreement between theory and observations. 

In connection with new observations further applica- 
tions of the evaporation model can be discussed (Liu, 
Meyer and Mcycr-Hofmoistcr, in preparation). Recently 
the concept of advection-dominatcd accretion was also ap- 
plied to the low- luminosity active galactic nuclei M81 and 
NGC 4579 (Quataert et al. 1999), NGC 4258 (Gammie et 
al. 1999) and low-luminosity elliptical galaxies (Di Mat- 
teo et al. 1999, 2000). It is not clear why for accretion 
rates which are different by about less than a factor of ten 
(scaled to Eddington accretion) the thin disks are trun- 
cated at > 10^ Schwarzschild radii or reach inward to « 
100 Tr-mS. Quataert et al. (1999) argued that the differ- 
ence might be caused by a different way of mass supply 
for the ADAF. 

9.4- Spectra from the standard thin disk plus corona 

One might ask whether spectra computed according to 
our modelling could be used to test our predictions. Such 
spectra would include radiation from an inner ADAF re- 
gion and an outer standard thin disk with the transition at 
a radius according to our M — r relation. The spectra de- 
termined according to the ADAF model (see e.g. Narayan 
et al. 1998) use the observational constraints for the inner 
edge of the cool disk and fit the observed spectra. If the 
derived rate M and the inner disk location agree with our 
M — r relation these already computed spectra are the 
spectra that our model will yield. 

10. Discussion of further aspects 

10.1. Magnetic fields 

The magnetic nature of friction in accretion disks is now 
generally accepted, supported by numerical experiments 



that show the generation and sustainment of magnetic 
fields in accretion disks of sufficient electrical conductiv- 
ity (for a recent study see Hawley 2000). This raises the 
question whether the neglect of magnetic fields in mod- 
elling the corona-disk interaction may not miss important 
eff'ects. 

Such effects arc twofold: (a) Where gyro-frequency is 
large compared to collision frequency electron thermal 
conductivity is only high along the magnetic field. This 
can significantly reduce effective thermal conductivity if 
the magnetic field is tangled. This can be taken into ac- 
count by reducing the conductivity coefficient Ko by the 
average of cos'^9 over horizontal surfaces, where 9 is the 
angle of the magnetic field with respect to the vertical, 
(b) Magnetic pressure contributes to the pressure sup- 
port. This can be taken into account by increasing the 
value 3?//i by a factor 1-|-^ in our equations, where the /3 
is the plasma parameter (which may also include turbu- 
lent pressiue). 

We have performed an order of magnitude approxima- 
tion to the individual terms in the global energy balance. 
As we will we show in separate work this allows to model 
our numerical results and suggests the following depen- 
dence on the parameters a, kq and /3. 
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where Tq is the saturation temperature reached at small 
radii and Tf is the subvirial temperature attained at large 
radii. The max;imum mass flow rate occurs where the two 
temperatures Tf and Tq become about equal. This yields 
the further proportionalities 
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This may serve as a guide on how the results depend on the 
choice of these parameter values. For values of /3 discussed 
in the literature see Sect. 8. 



10.2. The a- dependence 

For our numerical computations we took Q!=0.3. Values 
discussed in the literature include 0.1 to 0.3 from spectral 
fits of ADAF models to soft X-ray transients (Narayan 
et al. 1996), 0=0.25 obtained from applying detailed cal- 
culations to luminous black hole X-ray binaries (Esin et 
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al. 1997); and a=0.1 obtained from global magnetohy- 
drodynamical simulations of accretion tori (Hawley 2000). 
The above scaling relation shows how our results change 
if other a-values are chosen. 

10.3. Irradiation 

X-ray irradiation can cause a corona above the thin disk 
(e.g. London et al. 1981). This can interfere with the as- 
sumptions underlying the present model. Wc note that 
optically thick scattering layers in the inner corona can 
prevent such radiation to reach the disk surface (e.g. see 
Schandl & Meyer 1994). This can become a complex sit- 
uation which requires a detailed analysis to clarify under 
what circumstances this affects the applicability of our 
model. In general, models like ours should be good for 
black hole systems with ADAF type internal accretion. As- 
pects of formation of a disk corona have been investigated 
by de Kool & Wickramasinghc (1999). These calculations 
are more detailed with respect to ionization equilibrium 
than ours, but do not include thermal conduction, side- 
wise advection and wind loss in the energy balance. It is 
therefore difficult to compare with our results. 

11. Comparison with Honma's model 

Here we have considered the spatial co-existence of disk 
and corona and its gradual fading into a corona/ ADAF 
only region at its inner boundary. Honma (1996) in an 
apparently very different approach considered the effect 
of a turbulent diffusive heat flow outwards from a hot and 
mainly non-radiative advection-dominated inner region to 
an outercool accretion disk. In his simplified modelling 
both regions arc treated as one disk with radially varying 
thickness, and optical depth. He found a steep transition 
in the disk thickness and temperature at a radius that de- 
pends on the mass flow rate. In this radially thin transition 
layer most of the systems luminosity is radiated. 

Though at first sight very different, the two approaches 
have a basic physical similarity, though in an inverted ge- 
ometry: What is vertical in the disk-|-corona system is ra- 
dial in Honma's approach. Heat generated by friction is 
conducted outward (downward) from a hot inner (upper) 
region to a cool outer (lower) region where it is radiated 
away. In both cases this establishes the condition where 
the cool disk ends and the hot, ADAF type flow takes 
over. We note that Honma's diffusive heat conductivity 
bears resemblance to the thermal electron heat conduc- 
tivity. The rapid increase of the cooling with density leads 
to a similarly steep temperature profile in Honma's model 
as at the base of the corona in the disk-l-corona model. 

Wc have calculated the ratio of Honma's effective tur- 
bulent heat conductivity to that of electron thermal heat 
conductivity and obtain the value 1/6 for pressure and 
temperature at the peak of the coronal mass flow rate in 
Table 1. 



We think that the evaporation model is a more real- 
istic description of the transition from the cool disk to 
an corona/ ADAF only flow. In Sect. 2.5 we estimated the 
radial heat inflow and showed that it is a minor contribu- 
tion. Generally, the turbulent outward diffusion of specific 
energy reduces the inward advective transport (which is 
also diffusive) of energy. Thus the inward flow of energy 
is somewhat less than the inward flow of mass times the 
specific energy. Our approximate treatment of the inward 
flows of mass and energy in the derivations of Eqs. (39) 
and (46) in effect leads to such a distinction, for radial 
power laws as discussed by a factor of 2. 

In spite of the very different geometry and simplifi- 
cation Honma's model captures the same physical effect. 
His conclusion that his results shed light on the soft-hard 
transitions of black hole accretors is thus substantiated 
by the disk-|-corona model, though quantative differences 
exist. In Sect. 5 we have already remarked on where the 
"strong ADAF proposal" agrees with and where it de- 
viates from the results of the disk-|-corona model. This 
brings the "strong ADAF proposal", Honma's diffusive 
heat transport model, and the disk-|-corona model under 
a common physical aspect and makes a surprising similar- 
ity of some of their results results understandable. 

12. Conclusions 

Our model describes the interaction between the geomet- 
rically thin accretion disk and the corona above. We have 
derived simplified equations, which include the relevant 
processes determining the accretion via the cool disk to- 
gether with the hot coronal fiow which yields the matter 
supply for the ADAF towards the black hole. We deter- 
mine the structure of the corona which is in equilibrium 
with the cool disk. The equilibrium demands the evapora- 
tion of a certain amount of matter from the cool disk and 
we compute the evaporation rate as a function of the dis- 
tance r from the central black hole. This gives the rate of 
accretion of matter towards the inner region (part of the 
matter is lost by a wind). The evaporation efficiency be- 
comes higher with decreasing r, but reaches a maximum 
at some hundred Schwarzschild radii. For low accretion 
rates in the cool disk the evaporation uses up all matter 
at the distance where the evaporation rate is equal to the 
mass flow rate in the cool disk. The standard thin disk 
emits a near black body spectrum and, for high enough 
accretion rates the disk extends inward to the last stable 
orbit, the total radiation is soft. But if, for lower accretion 
rates the thin disk is truncated, the coronal flow/ADAF 
in the inner region provides a hard spectrum. 

These features of disk evaporation apply to galactic 
black hole X-ray binaries (X-ray transients, high-mass X- 
ray binaries) and also to AGN (truncated disks were found 
in low- luminosity AGN and elliptical galaxies). The mass 
accretion rates derived by spectral fits based on the ADAF 
model and the location of the inner disk (deduced from 
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observations or determined by the fitting) allow to check 
the predictions of the model for the location of transition 
to the coronal flow. These tests show good agreement for 
most cases of X-ray transients and the features connected 
with the maximum evaporation efficiency (see Sect. 9). We 
want to point out that our model gives qualitative results 
for relations between mass accretion rate, location of disk 
truncation and spectral transitions rather than exact val- 
ues. The division in an outer advection-dominated regime 
and an inner radiation-dominated regime of the coronal 
structure gives insight to basic features. Computed val- 
ues for the edge location might differ by a factor of two 
to three from the "true" values. An investigation of the 
detailed effect of simplifications in the equations as well 
as the effect of the a-value and of magnetic fields will be 
carried out later. 

A truncation of the thin disk farther in than 100 
Schwarzschild radii is in contrast to the predictions of 
our model. The reflecting component observed for Cyg 
X-1 seems to indicate such an inward extent of the thin 
disk. But, also for Cyg X-1, another investigation of 
observations leads to an inner edge location near 100 
Schwarzschild radii (see Sect. 9). This situation needs fur- 
ther clarification. But also from the theoretical side the 
possible existence of an "interior" thin disk below a corona 
in the innermost regions (closer in than the location of the 
evaporation efficiency maximum, compare Fig. 3) needs 
further investigation. 

For AGN it was found, based on ADAF spectral mod- 
elling that the truncation of the thin disk can be at very 
different distances from the supermassive black hole for 
mass accretion rates which differ much less (see Sect. 9). 
Our model only predicts the large distances. Also here ad- 
ditional investigations are needed to understand the phys- 
ical reason for such divergent results. 

Note added after submission of the manuscript: In 
a parallel investigation, very similar to our approach, 
Rozanska & Czerny ( ^stro-ph/0004158 , A&A in press) 
have also studied the coronal structure and the transition 
to an ADAF. They take into account the possibility of 
different temperatures of electrons and ions and include 
Compton cooling by the soft radiation from the under- 
lying disk. In our computed models for the parameters 
r/rg/^ > 25MEdd/M, that means log r/rg > 2.25 (com- 
pare Table 1), the energy equilibration time between elec- 
trons and ions is shorter than the heating time, thus es- 
tablishing temperature equilibrium between electrons and 
ions. For the same parameters also the rate of Compton 
cooling by the radiation from the underlying disk is small 
compared to the frictional heating rate and thus negligi- 
ble. 
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